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CZ2 ABSTRACT 

Aims. We present a detailed analysis of OGLE 2004-BLG-482, a relatively high-magnification single-lens microlensing event that exhibits 
^vq clear extended-source effects. These events are relatively rare, but they potentially contain unique information on the stellar atmosphere 
^ properties of their source star, as shown in this study. 

Methods. Our dense photometric coverage of the overall light curve and a proper microlensing modelling allow us to derive measurements of 
T-H the OGLE 2004-BLG-482 source star's linear limb-darkening coefficients in three bands, including standard Johnson-Cousins / and R, as well 

as in a broad clear filter. In particular, we discuss in detail the problems of multi-band and multi-site modelling on the expected precision of 

our results. We also obtained high-resolution UVES spectra as part of a ToO programme at ESO VLT, from which we derive the source star's 

precise fundamental parameters. 

y—{ Results. From the high-resolution UVES spectra, we find that OGLE 2004-BLG-482's source star is a red giant of MK type a bit later than 
M3, with Tgff = 3667 ± 150 K, log^ = 2. 1 ± 1 .0 and an assumed solar metallicity. This is confirmed by an OGLE calibrated colour-magnitude 
diagram. We then obtain from a detailed microlensing modelling of the light curve linear limb-darkening coefficients that we compare to 
^ model-atmosphere predictions available in the literature, and find a very good agreement for the / and R bands. In addition, we perform a 

• i-H similar analysis using an alternative description of limb darkening based on a principal component analysis of ATLAS limb-darkening profiles, 
and also find a very good agreement between measurements and model predictions. 

Key words, techniques: gravitational microlensing - techniques: high resolution spectra - techniques: high angular resolution - stars: atmo- 
sphere models - stars: limb darkening - stars: individual: OGLE 2004-BLG-482 



1. Introduction 

Photometric and spectroscopic observations of stars yield their 
specttal types and other information useful for studying their 
atmospheres, but much of the information on the structure of 
the atmosphere and related physical processes is lost in the 



* Partly based on observations made at ESO (073.D-0575A). 
** Royal Society University Research Fellow 



disc-integrated flux. Advanced models calculated for a broad 



range of stellar types (e.g. MARCS, Gustafsson et al. (20081; 



ATLAS, [Kimiczl ( [T992) ; |Plez et al.|(1992.) ) describe the corre- 
sponding physics at different optical depth, which can poten- 
tially result in observational signatures if the star's disc is spa- 
tially resolved. In particular, this information is present in the 
star's limb-darkening profile, which is the variation of intensity 
from the disc centre to the limb. Only a few observational meth- 



2 



M. Zub et al.: Limb-darkening measurements for a cool red giant 



ods such as stellar interferometry, analyses of eclipsing bina- 
ries, transiting extrasolar planets and gravitational microlens- 
ing are able to constrain stellar limb-darkening in suitable 
cases. Every single measurement thus provides an important 
opportunity for testing stellar atmosphere models. 

A Galactic gravitational microlensing event ( jPaczyhsE 



T986| occurs when a foreground massive object passes in the 



vicinity of the line-of-sight to a background star, resulting in a 
transient brightening of the source star (called magnification, 
or amplification). Microlenses can spatially resolve a source 
star thanks to caustic structures created by a lens. They are 
formed by a single point or by a set of closed curves, along 
which the point-source magnification is formally infinite, with 
a steep increase in magnification in their vicinity. In practice, 
this increase is so steep that the characteristic length scale of 
the differential magnification effect is of the order of a fraction 
of the source star's radius. Early works by e.g. Wittj( |1993] ) or 
Loeb & Sasselov ( 1995(1 have pointed out the sensitivity of mi- 



crolensing light curves to limb-darkening, with the aim to help 
remove microlensing model degeneracies. The specific use of 
microlensing as a tool to study stellar atmosphere was proposed 
later ( e.g. |Wlls-Gab aud||19^ |Sasselov|[l^ [Hendry et al.] 
19981 IGaudi & Gould 199"9ll, in particular to probe Galactic 



bulge red giant atmospheres ( Heyrovsky et al.||2000 i Indeed, 
for a given microlensing configuration, the spatial resolution 
increases with the source's physical diameter, so that giant stars 
are primary targets. 

Limb darkening measurements by microlensing were per- 
formed for a number of main-sequence and giant microlensed 
stars. Event MACHO 1998-SMC-l ( |Albrow et aT]|1999b| 
[Afonso et al.|[2000) l allowed for the first time such a mea- 
surement for a metal-poor A6 dwarf located in the Small 
Magellanic Cloud (SMC). Its stellar type was derived from a 
spectroscopic and photometric analysis in five filters; the lens 
was a binary star also located in the SMC. No real compari- 
son with atmosphere models could be provided since very little 
data existed for these metal-poor A stars. The first microlens- 
ing limb-darkening measurement for a solar-like star was re- 



ported by Abe et al. ( 2003| : the source was identified as an 
F8-G2 main-sequence turn-off star, involved in the very high- 
magnification microlensing event MOA 2002-BLG-33 caused 
by a binary microlens. A good agreement with limb-darkening 
coefficient predictions was obtained in the / band. A limb- 
darkening measurement for the late G / early K sub-giant was 



also performed by Albrow et al. (2001 1 with the binary-lens 
caustic-crossing event OGLE 1999-BLG-23. The stellar type 
of the source star was identified by comparing its position on 
two colour-magnitude diagrams obtained from two different 
telescopes, and deriving the star's effective temperature from 
colour calibration. Again, they found a good agreement with 
stellar models both for the / and R filters. 

Most of the limb-darkening measurements, however, were 
obtained on Galactic-bulge giant stars. The first case was re- 
ported by |Alcock et"aL] ( [T997] l for MACHO 95-30, which in- 
volved a very late M4 red giant source star (spectroscopic typ- 
ing). In this event theoretical limb-darkening coefficients were 
only used to improve the light-curve fit, but no limb-darkening 



argued that the intrinsic variability of the source star pre- 
cluded any useful limb-darkening analysis. Late M giants are 
of special interest because they give access to testing models 
at the lower end of the temperature range used to compute 
most of the synthetic model atmosphere grids. For the event 
MACHO 1997-BLG-28, 1 Albrow et aL] ( |1999c| ) derived / and V 
coefficients for a K2 giant (typing from spectroscopic observa- 
tions) crossing a caustic cusp, and found a good agreement with 
stellar models predictions. However, in such a complex event, 
many side effects could have affected the light curve, which 
somehow decrease the strength of the conclusions. Such a re- 
mark holds as well for MACHO 1997-BLG-41 (Albrow et al.l 



20001, which involved a late G5-8 giant crossing two disjoint 
caustics. 

Microlensing event EROS BLG-2000-5 provided the first 
very good opportunity to test at high precision the limb- 
darkening of a K3 giant (typing based on both photometry and 



high-resolution spectroscopy) in five filters (Fields et al. 2003 1. 
They concluded that their results in the V, I, and H filters were 
discrepant from atmosphere models, and furthermore argued 
that the discrepancy is unlikely to be due to microlensing light- 
curve modelling drawbacks, but could rather be explained by 
inadequate physics in the stellar models that may be not appli- 
cable for all surface gravities. A clear variation with time in the 
shape and equivalent width of the Ha line was also reported 
for the first time in this event ( |Afonso et al.|200r| [Castro et al. 
2001|). Limb-darkening was also detected in OGLE 2003-BLG- 



238 ( Jiang et al..2004) and OGLE 2004-BLG-262 (| Yoo et al. 
2004 1, which involved early Kl-2 giants, but no strong conclu- 
sions on limb darkening could be drawn from these events. 



From the binary-lens event OGLE 2002-BLG-069 (Cassan 



et al. 2004 Kubas et al. 2005 i, it was possible to obtain not only 
hmb-darkening measurements for a G5 bulge giant source star 
in the / and R bands, but also to directly test predictions from 
PHOENIX stellar model atmospheres by comparing the change 



of the Ha equivalent width during a caustic crossing (Cassan 



et al.[2004[[Thuri et al.[2006') using high-resolution UVESA^LT 



spectra. A discrepancy was found between model and observa- 
tions, which is most probably explained by the lack of a proper 
chromosphere implementation in the used stellar models. More 



measurement has been performed. Heyrovsky ( 2003 1 later 



recently, Cassan et al. ( 2006 1 performed limb-darkening mea- 
surements for the K3 giant source of OGLE 2004-BLG-254, 
and furthermore discussed an apparent systematic discrep- 
ancy between stellar model predictions and measurements that 
is observed for G-K bulge giants. However, in the case of 
OGLE 2004-BLG-254, it appeared that fitting all data sets to- 
gether or only a subset of them had an influence on the limb- 
darkening measurements ( [Heyrovs^ 2008| l, which remove the 
observed discrepancy. In order to quantify this effect, we pro- 
vide in this paper a detailed study on the impact of including 
data sets on the resulting limb-darkening measurements. 

We model and analyse OGLE 2004-BLG-482, a rela- 
tively high-magnification single-lens microlensing event that 
exhibits clear extended-source effects. The source-star funda- 
mental parameters and spectral typing were derived from a 
high-resolution spectrum obtained on VLT/UVES as part of a 
ToO programme. A good multi-site and multi-band coverage of 
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the light curve allows us to extract linear limb-darkening coef- 
ficients, which we compare to model-atmosphere predictions. 

The paper is organised as follows: in Sect. [2] we present 
the OGLE 2004-BLG-482 event, our photometric data and our 
data reduction procedures. We perform a detailed modelling of 
the light curve in Sect. [3] The fundamental properties of the 
target source star are derived in Sect. |4] Section 5 is dedicated 
to a detailed analysis of the measured linear limb-darkening 
coefficients and their comparison with model-atmosphere pre- 
dictions. Finally in Sect. [6] we perform a similar analysis using 
an alternative description of limb-darkening based on a prin- 
cipal component analysis of a set of ATLAS limb-darkening 
profiles. 



2. Photometric data 

2. 1. Event alert and follow-up observations 

The Galactic microlensing event OGLE 2004-BLG-482 
(a = 17h57m30.6s, 5 = -30°51'30"1 (J2000.0), or / = 
—0.3392°, b = —3.1968°) was discovered and publicly alerted 
on August 8, 2004 (MHJC [^ 3226 ) by the OGLE-Ii:0 Early 
Warning System ("EWS", Udalski|pOg3) l on the basis of ob- 
servations carried out in the / band with the 1.3 m Warsaw 
Telescope at Las Campanas Observatory (Chile). 



Table 1. Final selection of data sets, with the raw number 
of observational data (frames) and our final selection after the 
cleaning process. The last column lists the adopted error-bar 
rescaling factors. 



Telescope 


Filter 


Data 


(Frames) 




UTas (PLANET) 


I 


86 


(128) 


2.4 


Perth (PLANET) 


I 


13 


(15) 


3.8 


OGLE 


I 


44 


(68) 


2.4 


CTIO-Yale (jjPVN) 


I 


233 


(285) 


4.2 


Danish (PLANET) 


R 


51 


(67) 


3.2 


Auckland (^uFUN) 


(clear) 


266 


(334) 


2.4 


All data 




693 


(897) 





Following this alert, the PLANET collaboration (Probing 
Lensing ANomalies NETworIc) started its photometric follow- 
up observations on August 10, 2004 (MHJD ~ 3228), 
using a network of ground-based telescopes, including the 
Danish 1.54m (La Silla, Chile), Canopus Im (Hobart, 
Tasmania), and Perth/Lowell 0.6m (Bickley, Western 
Australia) telescopes. Data sets and quasi real-time fitted 
light curves were made publicly available onlin^ as part of a 
general data sharing policy. The event was also monitored by 
the /jFUN collaboratiorj^ which gathered data from six tele- 
scopes: the 1.3m and Yale 1.0m (Cerro Tololo Inter- American 
Observatory, Chile), the Palomar 1.5m telescope (Palomar 



' Modified Heliocentric Julian Date: MHJD = HJD - 2, 450, 000 
^ http://ogle.astrouw.edu.pl 
^ http://planet.iap.fr 



Observatory, USA), Wise Im (Mitzpe Ramon, Israel), and two 
New Zealand amateur telescopes at Auckland (0.35m) and 
Farm Cove (0.25m). 

On August 15, 2004 (MHJD ~ 3233), photometric data in- 
dicated a deviation from a normal point-source point-lens light 
curve. A public alert was issued on August 16, 2004 16:05 
UT, pointing towards a high peak magnification event, possi- 
bly featuring strong extended-source size effects. In the follow- 
ing hours, on August 17, 2004 a target of opportunity (ToO) 
was activated on the UVES spectrograph at ESO VLT in or- 
der to monitor the event peak magnification region where spec- 
troscopic effects are expected. Thanks to an almost real-time 
modelling operated in parallel, the crossing time of the source 
disc by the lens was estimated to be around 2.4 days. The peak 
of the light curve was reached on August 18, 2004 18:32 UT 
at almost three magnitudes above the baseline, corresponding 
to a minimum (i.e. with null blending) peak magnification of 
A - 15. 



2.2. Data reduction and error bars 

The OGLE data were reduced with their own pipeline, while 
PLANET and /jFUN data were reduced with various versions 
of the PLANET pipeline (pySIS; [Aibrow et al.|2009"l l. All these 
reductions are based on the image-subtraction method ( |Alard| 
& Lupton|1998 Bramich|2008) . A preliminary image-quaUty 
inspection helped to remove images with a significant gradient 
across the field, owing to strong background moonlight. Under- 
exposed images were also removed in this process. We paid 
particular attention to the quality of data taken at La Silla at the 
time of peak magnification, because of unfavourable weather 
conditions at that site. We could however keep a few trusted 
data points. 
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After the data reduction process, we set for each PLANET 
and /jFUN telescope a range of seeing and background within 
which the homogeneity of the data sets is ensured. For the Yale 
telescope, we excluded data with seeing outside the range 1.8- 
3.2". In the case of UTas data, we applied an upper limit on 
the seeing of 3.0", and for the Perth, Danish, and Auckland 
telescopes, 3.3". Our final data set is presented in Table[T]and 
displayed in Fig. [T]( Auckland telescope had no filter at the time 
of the observation, so the filter is referred to as clear). 

It is known that the error bars we obtain from the data re- 
duction are usually underestimated, and are not homogeneous 
from one data set to another. To avoid this problem, we rescaled 
the eiTor bars, so that from the best model one has X^/N ~ 1 
for each data set fitted alone, with N the coiTesponding number 
of data points. Moreover, it happens that some of the original 
eiTor bars o are unrealistically small; to prevent this, we added 
in quadrature an additional term to the rescaled eiTor bars, so 
that o'^ = (fcaO)^ + (4 X 10""*)^ magnitudes. When the origi- 
nal error bars are small (UTas and Danish), the constant term 
dominates the error bars' size. The values of Ica are given in 
Table [1] 
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3230 3235 3240 3245 3250 

t = HJD-2,450,000 



Fig. 1. Light curve of OGLE 2004-BLG-482, with data from PLANET (Danish, UTas and Perth), OGLE and ^FUN (CTIO-Yale 
and Auckland) collaborations. The two gray solid lines in the upper panel draw the best-fit model for the / and R filters with 
linear limb-darkening parameters given in Table |2] The two dotted curves correspond to the two extreme cases, F = (uniformly 
bright source, lower dotted curve) and F = 1 (upper dotted curve). The two pairs of vertical dashed lines marked m = p* and 
u = 0.5 p* indicate when the lens is located at the limb of the source and half way from its centre to the limb. All the curves 
intersect at m = 0.77 p*, also marked by a vertical dashed line. The fit residuals in magnitudes are displayed in the lower panel. 



3. Light-curve modelling 

3.1. Linear limb-darkening formalism 

Limb-darkening profiles of stars can be described analytically 
at different levels of approximation, in particular by a sum con- 
taining powers of /J = cos a, where a is the angle of a given 
emerging light ray with respect to the normal of the stellar sur- 
face (e.g. |Claret 2000 1. In the first degree of approximation, 
called the linear limb-darkening (hereafter, LLD) law, the star 
brightness profile can be written as 



I{r) = l-a (^1- - 



(1) 



where r = ^/F— /j-" 
disc from where the light is emitted, and a is the linear 
limb-darlcening coefficient (hereafter LLDC). In this work, we 
will concentrate on measuring LLDCs. Firstly, because in mi- 
crolensing events higher order coefficients have a very small 
impact, e.g. |Dominik| pO(M| finds that for a caustic crossing, 
the effect of the change of the LLDC on light curve is ^ 25 
times greater than the square-root coefficient. Secondly, be- 
cause a strong correlation exists between the coefficients, it is 
impossible to precisely measure the LLDC when a further coef- 
ficient is taken into account (I Kubas et al.|2005| l. Lastly, because 
LLDC are widely used and are available in catalogues; it is an 
important aspect for our goal to compare our results with the 
existing literature. 



For our modelling purpose, a more convenient way to 
rewrite the LLD law is to have a formula that conserves the 
total source flux for all LLDC values. With this requirement, 
the LLD law equivalent to Eq. ^ (Albrow et al.||I999a| l but 
normalised to unit flux can be written as 



n 



1-F 



where F is the LLDC modelling parameter, with 
3F 



2 + F 



(2) 



(3) 



With this formalism, it is interesting to notice that all limb- 
darkening profile s intersect at a common fractional radius 
|Heyrovsky|2003|), rm = V5/3 ~ 0.75. 



3.2. Single-lens, extended-source models 

In its motion relative to the lens, the source centre approaches 
the lens at a minimal distance mq in units of the angular Einstein 

ring radius = ^4GMc-^{Dl^ -D^^) < Einstein|l936 wifli 
Ds, Dl the distances from the source and the lens to the 
observer, M the lens mass), which can be smaller than the 
source radius p* expressed in the same units. Because high- 
magnification events involve low values of the impact parame- 
ter Mo, they are likely to be affected by extended-source effects 
in particular if the source star is a giant. Although this happens 
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fairly rarely in practice (a couple of cases amongst the ^ 700 
microlensing events observed every year), this is the case for 
OGLE 2004-BLG-482. 

The point-source magnification at the exact location of the 
lens is formally infinite, following the well-known formula 
( |Paczyriski|1986 i 

+2 

'^p.splW^^^^r, (4) 
M VM-^ +4 

where u is the distance from the lens to a given point on the 
source in units of Og. Consequently, the flux originating from 
regions of the source in the immediate neighbourhood of the 
lens (typically a fraction of the source radius) is preferentially 
amplified. The relative motion of the source and lens then re- 
sults in a time-dependent probing of the stellar atmosphere at 
different fractional radius, corresponding to different optical 
depths. 

Single-lens light curves affected by extended-source effects 
display a characteristic flattening at their peak. For a uniformly 
bright extended source, [Witt & Mao| ( |1994[ ) derived an exact 
analytic formula for the magnification, which involves ellip- 
tic integrals. But there is no similar formula to describe limb- 
darkened sources, and calculating the exact magnification re- 
quires numerical integration. One way is to decompose the 
source into small rings of uniform intensity. Another approach 
by Heyrovsky ( 2003[ l is to perform the angular integration over 
the stellar disc analytically and only the radial integration nu- 
merically, for arbitrary sources. 

If some conditions are fulfilled, it is also possible to use 
approximate formulae, which have the advantage to allow us 
a very fast computation. Considering that in Eq. Apspi ~ 
1/m when m <C 1, Yoo et al. (20041 find that the magnification 



A\\i for an extended source with a linear limb-darkening profile 
with coefficient Y can be expressed as 

Aiid(M,p*) = [Bo(z)-rBi (z)] Apspi(M) , 
z = m/p* , 



4z 

71 



arcsinmin ( 1, - ) ,z 



(5) 



3z 



71 1 



^ ^ \J ^7} — 2zrcos(|) 



r.<l), 



where E is the incomplete elliptic integral of the second kind 



following the notation of Gradshteyn & Ryzhik (1965i. The 
integral B\ can be efficiently evaluated and tabulated for z, as 
can Bo- This approximation is valid as far as p^/8 <C 1 and 
Mo ^ L Because these relations hold for OGLE 2004-BLG- 
482 we choose this formalism (although close to the limit case 
of application, since the maximum error for a uniform source 
here is of the order of 0.2%, but is still much lower than the 
photometric errors). 

The complete model then involves four parameters: the 
source radius p*, as well as mq, fo and Je, which define the 
rectilinear motion of the source with respect to the lens, so 
that the lens-source separation u satisfies u^{t) = Mq + (f — 
tof' /t^. Moreover, one has to take into account for each tele- 
scope "/" two more parameters, the baseline magnitude M]^ — 



—2.5 log(fg +F^) and the blending factory' =F^/F^. Here, 
and are the source and the blend flux, the latter referring to 
any un-magnified flux entering the photometric aperture, from 
the lens itself and e.g. background stars. They are related to the 
time-dependent magnification A\\i by F'{t) — Am {t)F^+F^. 



3.3. Fitting procedure 

To fit our data sets, we use two minimisation schemes: Powell's 
method and a Markov-Chain Monte-Carlo algorithm, from 



which we also obtain the model parameter error bars (Kains 
|et al.||2009 Cassan et al.||2010 i. As stated before, it is impos- 
sible to define a proper number of degrees of freedom. Indeed, 
the parameters uq, fo, fe and p* are common to all data sets, 
whereas Mj^ and g' are associated to the data set and the 
LLDCs may be chosen to be common per observing filter or per 
individual telescope. This explains the choice of instead of 



d.o.f to rescale the error bars in Sect. 2.2 The first requirement 
to get precise measurements of limb-darkening coefficients is 
to get an overall well-covered light curve. This allows us to se- 
cure good measurements of the basic parameters uq, fo, fe and 
Ph., as well as and g'. The region of the light curve with 
a noteworthy sensitivity to Umb-darkening is, however, mainly 
limited to when the lens is inside the source-star disc, and drops 
to a few percent outside ( Heyrovsky et al.|2000 i. We now dis- 
cuss this aspect in greater detail. 

While all limb-darkening profiles intersect at the same frac- 
tional radius r\\i ~ 0.75 as seen in Sect. |3.1[ the correspond- 
ing magnification light curves intersect at around Miy ~ 0.77 p^, 
(with u the lens-source centre distance). This special point is 
marked by a vertical dashed line in Fig. [T] in which we have 
also indicated two other interesting positions of the lens: at the 
limb of the source (m = p*) and at half-way from its centre to 
its limb (m ~ 0.5 p*). The two dotted magnification curves of 
the figure show the two extreme cases of LLDC, F = (no 
limb-darkening) and F = L From this we can distinguish three 
main regions: < m/p* < 0.5, where the limb-darkening sen- 
sitivity is high, up to ^ 16% ; 0.5 < m/p* < 0.77, where the 
sensitivity decreases outward to 0, and 0.77 < m/p* < 1 where 
the sensitivity increases outward and peaks at the limb at ~ 8% 
( |Heyrovsky|2003| l. Based on this argument and from our data 
coverage of OGLE 2004-BLG-482 shown in Fig. [T] it is clear 
that we can expect LLDC measurements from UTas /-band, 
Danish /?-band and Auckland's clear-filter 

The best-fit parameters and their error bars are given in 
Table |2] for different combinations of data sets. We comment 
on the results in detail in Sects [5] and [6] In Fig. [T] we plot the 
combined fit including all telescopes and using one coefficient 
per band. 



3.4. Estimates of tlie lens properties 

Although the properties of the lens are not of primary interest 
here, we can still provide an estimate of the lens' mass and dis- 
tance. However, these quantities cannot be measured here, be- 
cause an additional observable, such as the parallax, is needed 
to remove a degeneracy between these two parameters. Here, 
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parallax effects are not visible because the time scale of the 
event is very short, ~ 10 days ^ 1 year 

From our modelling and our estimate of the source radius 
and distance (see Sect. |4.1| l, we derive the Einstein radius to be 
around 9e = 9*/p* ~ 0.4 mas, which leads to a relative proper 
motion of /j = Gn/fE — 16 mas/yr. This high proper motion al- 
most certainly means that the lens is located in the disc (or 
possibly in the halo). Moreover, with such a high motion, there 
is a good chance that the lens can be clearly visible (away from 
the source) in a few years by adaptive optics observations. 
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16 



o 



20 




Fig. 2. OGLE-m BLG182.8 field calibrated / vs. {V - 1) 
colour-magnitude diagram, comprising stars within a radius of 
2.16' centred on our target OGLE 2004-BLG-482 (red open 
triangle). The red circle indicates the mean position of the 
RCG centre, and the cross the width of the two-dimensional 
Gaussian distribution. 



4. Source star properties 

4. 1. OGLE calibrated colour-magnitude diagram 

The microlensing event OGLE 2004-BLG-482 occurred in 
OGLE-III BLG182.8 field, and was also observed during the 
second phase of OGLE in field BUL_SC23. From the calibrated 
photometry in / and V filters of the OGLE-III BLG182.8 field, 
we extract an / vs. {V — I) colour-magnitude diagram (CMD) 
by selecting stars surrounding our target within a circle of ra- 
dius 2.16' (~ 9000 stars), as shown in Fig. [2] This choice en- 
sures that we have enough stars to construct the CMD while 
keeping a reasonably homogeneous extinction across the se- 
lected region. 

Our target is indicated as the red open triangle and has 
a calibrated magnitude and colour of / = 14.41 ± 0.03 and 
(y -/) = 3.93 ±0.04. From the analysis of the OGLE im- 
ages, we conclude that this bright target is not blended by a 
neighbouring star in /. We also checked in the OGLE-III photo- 
metric catalogue, which has a better resolution than OGLE -II, 
that there is no blended star within 1 arcsec bright enough in 



V to contaminate our measurements. Finally, our model finds 
a blending ratio close to zero, justifying our assumption that 
the measured magnitude and colour of the target can safely be 
assigned to the source star 

In order to correct these measurements for extinction, we 
can either use the reddening maps of 'Sumi| ( |2004| l based on 
OGLE-II photometry, or use the red clump giant (RCG here- 
after) assuming that our source suffers the same amount of 



extinction. At the position of the source, Sumi (2004j| mea- 
sures an extinction of ^(y — /) = 1.405 ±0.027. This is de- 
rived assuming a clump colour of 1.028 and a ratio of to- 
tal to selective extinction Ri = 0.964, giving an absorption of 
Ai = /?i X £ (y - /) = 1 .36 ± 0.06. 

The RCG central position is marked in Fig. |2] as a red 
circle with error bars. To determine its mean magnitude and 
colour, we fit a two-dimensional Gaussian around its position 

400 stars), from which we derive /rcg = 15.88 ±0.01 and 
(y — /)rcg = 2.263 ± 0.004. Using the same clump colour as 
|Sumi| ( 255^ , we get an extinction of E{V -I) = 1.235. Given 
the uncertainty of the intrinsic clump colour, due to variations 
with age and metallicity ^Salaris & Girardi|2002 1, this estimate 
agrees with the previous one. We therefore adopt as the dered- 
dened magnitudes and colour of our target Iq = 13.05 ±0.07, 
(y - /)o = 2.53 ± 0.05 and Vq = 15.58 ± 0.09. 

In principle the observed position of the clump could be 
used to measure its distance. In practice, the absolute magni- 
tude of the clump depends on age and metallicity, and cor- 
rections introduce an uncertainty as for its intrinsic colour. 
Moreover, the reddening corrections are not accurately deter- 
mined. We therefore prefer fixing the distance and check that 
the observed clump position is compatible with this choice. We 
start by assuming a Galactic Centre distance of 8.33 ±0.35 
kpc from Gillessen et al.| ( |2009| l, giving a distance modulus 
of 14.60 ± 0.09. Then, we assume that Baade's Window at 
(/ = 1.00° b = —3.88°) is at about the same distance as the 
Galactic Centre, according tolPaczynski & Stanek ( 1998 1. This 
fixes the distance to OGLE-II field BUL_SC45, which contains 
Baade's Window. Finally, we use the relative distance of field 



BUL_SC23 with respect to BUL_SC45 as given by |Rattenbury 
et al. ( 2007| l (their Table 1), which amounts to 0.13 mag, to 



get a distance modulus of 14.73 ± 0.15 to the clump in the di- 
rection of our target. The corresponding absolute magnitude of 
the clump is then 15.88— 1.36— 14.73 = —0.21, in good agree- 
ment with the most recent value determined by |Groenewegen| 
(j2008j) for the local red clump based on revised Hipparcos data, 
namely M/,rcg = -0.22 ± 0.03. 

If our source is at the same distance as the clump in its di- 
rection and suffers the same amount of extinction, its expected 
dereddened magnitude is 14.41 - 15.88 + 14.51 13.04, in 
good agreement with the previous estimate based on Sumi's 
reddening law. The agreement in colour is not as good, at 
3.93-2.263 + 1.028 = 2.70. 

We then fit calibrated isochrones from |Bonatto et al.]p004| l 
to 2MASS data in our field, to derive the following near- 
infrared extinctions: Aj = 0.52 ±0.10, Ah = 0.36 ±0.11 and 
Ak^ — 0.20 ± 0.02. From this and the magnitudes listed in 
the 2MASS PSC for our target (2MASS 17573061-3051305), 
we get 70 = 11.55 ±0.10, Hq = 10.68 ±0.11 and K^fi = 
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Fig. 3. Observed (black line) and best-fit (blue) template spectrum of OGLE 2004-BLG-482. The region around the TiO 7100 
shows the agreement of the observed and synthetic spectra. The two orange curves are plotted at ±100 K. The regions excluded 
from the fitting process are shaded in light grey, while the remainder of the spectrum is shaded in dark grey. 



10.42 ±0.04, and the coiTesponding colours {J-H)o = 0.87 ± 
0.16, {H-Ks)q = 0.26±0.12 and {J - Ks)o = 1.13±0.11. 
Converting to Bessell & Brett near-infrared photometric system 
( [Bessell & Brett|1988| gives {J-Ks)q = 1.17 and /To = 10.46, 
corresponding to an M4 giant (their Table 3), which have mean 
colours of (y — /)o — 2.55 and {V ~K)() = 5.10, in good agree- 
ment with our observed values (V — /)o = 2.53 and {V —K)q = 
15.58-10.46 = 5.12. 

This allows us to estimate the source radius using the sur- 
face brightness relation: log9, ± go/5 = (0.037 ± 0.0 07) x 
{V - K)q + (0.610 ± 0.028) from 'Groenewegenl (|2004| cali- 



brated on 40 M giants, where 9* is the source angular diameter 
in mas. We find an angular diameter of 9* = 51 ± 3 /jas, which 
at the adopted source distance of li = 8.8 ±0.6 kpc corresponds 
to a physical source radius of R^, ~ 48 ±4 

In the next section, we perform the analysis of the 
VLT/UVES high-resolution spectra that we obtained on this 



event, in order to derive more accurately the spectral type and 
to determine the fundamental parameters of the source star. 



4.2. VLT/UVES spectroscopy 

We have obtained for OGLE 2004-BLG-482 high-resolution 
spectra (R - 40000) on VLT/UVES, as part of a ToO activated 
shortly after the peak of the light curve was passed. The data 
were reduced in a standard way using version 2.1 of the UVES 
context of the MIDAS data reduction software. 

The spectrum is dominated by broad absorption bands 
from molecules. The shape and depth of molecular absorption 
bands, particularly TiO, are very sensitive to the stellar effec- 
tive temperature T^ff, and to a lesser degree also to the sur- 
face gravity \ogg. We estimated the atmospheric parameters 
of OGLE 2004-BLG-482 by comparing the observed spectrum 
with a grid of pre-calculated synthetic template spectra. 
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The grid of synthetic template spectra, calculated by Plez 
(priv. comm.), is based on synthetic spectra calculated from 
MARCS spherical model atmospheres with ID emergent spec- 
tra and LTE radiative transfer ( .Gustafsson et al.|[2008l [2D03] 
[1975) |Plez et al.|2003l [T992| , a nd includes the latest available 
atomic and molecular hne data d Gustafsson et al.|2008[ |Kupka| 
et al.|1999||Plez|1998] l. Synthetic template spectra for M giants 



calculated with the MARCS model atmospheres have a good 
record for determining stellar parameters in M supergiants (e.g. 
Levesque et al.|2005 2007 Massey et al.|2008j ) and were ex- 



tensively used to calibrate M giant photometry ( [Bessell et al.| 
[1998). 

The grid used in our analysis covers an effective tempera- 
ture range of 3000K < T^ff < 4000K, with steps of 100 K, and 
a surface gravity range of 0.0 < \ogg < 3.0, with steps of 0.5. 
This grid was calculated for giants with solar abundances and 
no carbon enrichment. Since our grid does not cover a range 
of metallicities, we therefore have no leverage on this parame- 
ter. We also prepared routines to calculate linear interpolations 
between the spectra in our grid for any given value of Teff and 
log§- 

We then compared the observed spectrum of OGLE 2004- 
BLG-482 with template spectra across the available range of 
Tgff and \ogg and determined the goodness-of-fit using the 
diagnostic. In calculating x^, we used the entire observed spec- 
trum, from approximately 4800 to 10000 A, only excluding 
three regions that are strongly affected by telluric absorption 
(7580-7850, 9070-9120A and 9300-9800 A). However, since 
no continuum is present in our spectrum, and we also do not 
know the absolute stellar flux, we renormalised the synthetic 
spectrum using a one-dimensional polynomial function prior to 
calculating x^. This renormalisation does not affect the shape 
of the broad molecular bands that are important for determining 
Jeff and log^. 

We refined the best values of Teff and log^ using parabolic 
minimisation between the grid points that yielded the lowest 
value of the diagnostic. In Fig. [3] we illustrate the agree- 
ment between the observed and best-fit template spectrum, in- 
cluding estimated parameter uncertainties, around the highly 
temperature-sensitive TiO band near 7 100 A. We find that the 
parameters that best fit our observed spectrum are T^-ff = 3667 ± 
150K and \ogg = 2.1 ± 1.0, assuming solar abundances. The 
quoted error bars are dominated by systematic uncertainties in 
the synthetic spectra and data reduction procedures used, such 
as flux calibration. Our uncertainties are further increased be- 
cause our grid of template spectra was calculated for only one 
metallicity. The range of effective temperatures we find is com- 
patible with a star of MK spectral type between Ml and M5, 
with the best-fit value giving a red giant star a bit later than M3 
( |Houdashelt et al.|2000[|Strassmeier & Schordan|2000| l. 

The large error bar on the surface gravity confirms that 
our spectrum has little to offer in gravity-sensitive diagnos- 
tics. However, we can obtain independent constraints on \ogg: 
given that the mass of an M giant of 1 or 10 Gyr is smaller than 
2.3 and IM0 respectively, using logg — logg^ +logM — 2 x 
logR^,, we find the corresponding upper limits of the surface 
gravity: log^ — 1.5 ±0.2 and log^ =1.1 ±0.2 respectively, 
taking into account the uncertainty of the source radius. This 



agrees with our spectroscopic analysis, although favouring the 
lower boundary. 

4.3. Conclusion on the source MK type and parameters 

We finally find a good agreement between our photometric and 
spectroscopic study, with a source star of MK spectral type 
a bit later than M3. We therefore adopt the fundamental pa- 
rameters from the spectroscopic analysis (Teff = 3667 ± 150 K, 
logg = 2.1 ± 1.0, solar metallicity) to make our selection of 
atmosphere models used to compare our limb-darkening mea- 
surements to model predictions, as discussed in the next sec- 
tion. 

5. Linear limb-darkening discussion 



As discussed in Sect. |3.3[ three data sets have some sensi- 
tivity to limb-darkening: UTas (/-band), Danish (/?-band) and 
Auckland (clear filter). The first question we address now is 
how the individual linear limb-darkening coefficients (LLDC) 
are affected by including or removing some of our data sets. 
Indeed, our first step was to model every data set indepen- 
dently, and step by step to include other telescopes. We first 
noticed that there was a change in the LLDC values that de- 
pends on the added data sets. We therefore performed a detailed 
analysis to understand what could cause such variations, and 
to identify combinations of data sets that lead to correct LLD 
measurements. The results we are commenting are presented 
in Fig. |4] the three columns correspond to UTas, Danish, and 
Auckland, respectively, and the individual panels display the 
LLDC measurements for various combinations of data sets; the 
corresponding model parameters are given in Table [2] In the 
figure and table, the letters A, C, D, O, P, U refer to the tele- 
scopes Auckland, CTIO Yale, Danish, OGLE, Perth, and UTas 
respectively. 

UTas (U) clearly provides the best data set for LLDC mea- 
surements, since the data sample the whole LLD-sensitive re- 
gion at the peak of the light curve, as well as its wings and 
baseline. On the other hand, modelling Danish (D) alone pro- 
vides a very unrealistic fit, with large error bars and very ir- 
regular MCMC correlations. To explain this, we recall that as 



mentioned in Sect. 2.2 the peak of the light curve was observed 



under bad weather conditions in La Silla, in particular the two 
consecutive data points around t = 3235.5. Moreover, the data 
coverage is not optimal, because there are only two epochs that 
cover the LLD-sensitive part of the light curve. As a result, this 
poor coverage combined with some uncertainty in the data lead 
to a model that apparently looks better in terms of chi-squared, 
but cannot be trusted. The last telescope with data sensitive to 
limb darkening is Auckland (A). We can fit the corresponding 
data alone and obtain a reasonable fit, but we obtain large er- 
ror bars because the photometric accuracy of the data is several 
times lower than for UTas, and furthermore, the data taken dur- 
ing the source crossing are all located close to the limb, in the 
region of less sensitivity to limb darkening. We note that the 
LLDC we obtain is higher than UTas's, which is expected, be- 
cause Auckland's clear filter is known to peak between red and 
infrared and LLDCs usually decrease towards the infrared. 
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Table 2. Model parameters and error bars for different relevant combinations of data sets. The measured linear limb-darkening 
coefficients are indicated in bold face. The data sets are referred to by letters, following the convention indicated in the first line 
of the table. The uncertainties on the parameters are indicated in parenthesis and apply to the last significant digit. Models for 
which no stable fit or very unreaUstic results are obtained are marked with the symbol " [★]" following the measured value. 



Parameters 


UTas (U) 


Danish (D) 


Auckland (A) OGLE (O) Perth (P) 


CTIO Yale (C) 


Independent fits for V, D and A 


to (days) 

«0 

tE (days) 

p* 

a 

Mb 
g 


3235.78(4±1) 

0.010(8±4) 

8.(9±1) 

0.14(0±2) 

0.677 ±0.013 

11.5 

7.0 

82.5 


3235.78(3 ±3) 

0.0(2±1) 

9.(6±4) 

0.1(3±1) 

0.67 ± 0.22 [★] 

11.4 

1.4 

43.2 


3235.76(8 ±4) 

0.0(0±1) 

9.(3±3) 

0.14(0 ±7) 

0.76 ±0.13 

13.5 

4.5 

234.9 




Combined fit including U+D 


to (days) 
«0 

tE (days) 

P* 

a 

Mb 
g 


3235.784(5 ±8) 

0.00(9 ±2) 

9.1(5±9) 

0.13(7±1) 

0.674 ±0.012 

11.5 

7.2 

85.1 


0.837 ±0.018 

11.4 

1.3 

57.8 






Combined fit including U+A 


to (days) 

UQ 

tE (days) 

p* 

a 

Mb 
g 

x' 


3235.780(8 ±8) 

0.00(0 ±3) 

9.(1±1) 

0.13(8±2) 

0.714 ±0.013 

11.5 

7.2 

101.9 




0.660 ±0.023 
13.5 

4.5 
286.3 




Combined fit including D+A 


to (days) 

Uq 

tE (days) 

p* 

a 

Mb 
g 


3235.77(5 ±3) 
0.00(0 ±7) 
9.(7±2) 
0.13(4±8) 


I. ± 0.23 M 

II. 4 
1.4 
50.6 


0.93 ± 0.29 [*] 

13.5 

4.8 

241.5 




Combined fit including U+D+A 


to (days) 

UQ 

tE (days) 

pt 

a 

Mb 
g 

x' 


3235.781(4±9) 

0.00(0 ±4) 

9.2(9 ±6) 

0.13(6±1) 

0.713 ±0.012 

11.5 

7.3 

102.7 


0.881 ±0.010 

11.4 

1.3 

58.1 


0.660 ±0.011 

13.5 

4.6 

287.8 




Combined fit including all telescopes (one LLDC per band) 


to (days) 

"0 

tE (days) 

p* 

a 

Mb 

g 


3235.781 (6 ±7) 

0.00(0 ±2) 

9.6(1 ±2) 

0.130(9±5) 

0.714 ±0.010 

11.5 

7.6 

122.7 


0.884 ±0.021 

11.4 

1.4 

51.0 


0.652 ±0.016 0.714 ±0.010 0.714 ±0.010 
13.5 14.1 12.7 
4.8 0.0 0.7 
286.6 42.6 14.3 


0.714 ±0.010 
14.0 
-0.8 
239.7 
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Fig. 4. Graphical representation of the linear limb-darkening measurements (crosses) for the three data sets with sensitivity to 
limb-darkening: UTas in the /-band, Danish in the R-hwA, and Auckland in a clear filter. The open hexagons and the filled 



diamonds are the predictions from Claret ( 200()|) and Heyrovsky (2007 ) linear limb-darkening (LLD) coefficients. The fitting of 
the light curve is performed for different combinations of telescopes (same letter conventions as for Table |2]i, and the results are 
discussed in Sect.|5] The adopted measurements are those marked with black squares in the upper right of the panels. 



Starting from these models, we include different combina- 
tions of other data sets. If we base our analysis on the LLDC 
measurement from our best data set, UTas, then we find two 
distinct behaviours: either the UTas's LLDC is not displaced 
from the individual fit (a ~ 0.67, e.g. U or U-i-D) or is slightly 
modified (a ~ 0.71, e.g. U-i-A or all telescopes). Interestingly, 
when combining U and D data, the fit is stabilised for D. This 
is because the common fitting parameters (Mo,fo,fE, p*) are bet- 
ter constrained. However, combining A with U data modifies 
the LLDCs compared to A and U modelled alone. This is also 
true when combining U-i-D-i-A or all telescopes. That Auckland 
modifies UTas's LLDC (our best-suited data set for LLDC 
measurements) when the two data sets are combined, lead one 
to be careful about the interpretation of Auckland's LLDC (be- 
sides the large error bar on the LLDC). 



From this, we conclude that a precise measuring of LLDC 
requires a very careful study: first, one has to identify the data 
sets that can potentially provide a limb-darkening measurement 
with enough sensitivity, based on the light-curve sampling as 
discussed in Sect. |3.3| Then, one has to check whether the in- 
clusion of additional data sets affects the results. Indeed, as we 
have shown for this microlensing event, adding more data sets 
to the light curve modelling can lead to two opposite effects: 
either the new data stabilise the fit and help obtain LLDCs 
for more data sets, or they perturb the LLDC measurements. 
The latter may happen if unknown systematic effects are af- 
fecting the data. For OGLE 2004-BLG-482, the most reUable 
LLDCs for UTas and Danish are obtained when these data sets 
are combined in the fit (U-i-D). No definitive conclusion can be 
safely drawn for Auckland LLDC, although its best estimation 
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is likely to be obtained using A data alone. The relevant mea- 
surements we discuss below are marked in Fig. [4] with a black 
square in the upper right of the corresponding panels. When 



the fit is performed using the formula of Heyrovsky ( 2003 i, 
we obtain similar results for the combinations Uh-D and A: 



flu = 0.655 



hO.OlO 
-0.016' 



flD = 0.825 



h0.023 
-0.022 



and flA 



. 7^, +0.083 
^- '-^^-0.096- 



In order to compare our measurements to linear limb- 
darkening predictions from atmosphere models, we use two 
sets of LLDCs computed from Kurucz's ATLAS models, 
for which hydrostatic equilibrium and LTE were assumed 
[T9921 [1994}. The first set of LLDCs is taken 
using the VizieR database, for the whole 



Kurucz 



Claret 



2000) 1, 



(e.g. 
from 

available range of temperatures and log^ compatible with 
OGLE 2004-BLG-482's source star fundamental parameters 
(Sect. |4|i; we assume a solar metallicity to be consistent with 
our spectral analysis. The corresponding LLDCs are plotted in 
Fig.|4]as thin, open hexagons. We find twelve models that cor- 
respond to our requirements: two different temperatures (3500 
and 3750 K), three \ogg (1 .0, 1 .5 and 2.0, plotted from smaller 
to bigger symbols) and for each configuration two microtur- 
bulent velocities (1 and 2 km/s). The second set of LLDCs is 
plotted as filled diamonds, and correspond to coefficients com- 



only Heyrovsky (2007) predictions are available; but within the 
large error bars commented on previously, the data are fairly 
compatible with the predictions. 

6. PCA-based limb-darkening coefficients 

Although stellar limb darkening is usually modelled by ana- 
lytical laws, another option is to construct new bases of func- 
tions directly from model-atmosphere limb-darkening profiles. 
In this section, we use a limb-darkening basis numerically con- 
structed by principal component analysis (PCA and PCA LD 
in the following) for a set of given model atmosphere limb- 
darkening profiles, following [Heyrovs ky (2003). 

In this approach, the stellar intensity profile is expressed as 



/(r)=^a,/,(r), 



(6) 



in place of Eq. ([T]), where the /,(r) are the PCA basis func- 
tions and a,- the corresponding coefficients. For our analysis 
we used a very general PCA basis constructed from the BVRI 
profiles of the full Kurucz (1992) ATLAS model-atmosphere 
grid (see Heyrovsky 2008 for details). The resulting three first 
basis functions, computed for a set of atmosphere models that 
match the stellar parameters of OGLE 2004-BLG-482's source 
star, are displayed in Fig.|5] 

In the simplest case of a 2-term PCA LD law (the analogue 
of the analytical linear limb-darkening law, LED), the relevant 
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puted using the interpolation method advocated by Heyrovsky] 
( 2007| l. These are computed for the same stellar parameters as 
before. 

Evidently our LLDC measurements agree very well with 
the predictions from atmosphere models. For UTas /, our mea- 
surement is compatible with both the predictions from |Claret| 
( pOOO) ! and |Heyrovsky| (g007|. For the Danish R filter, the 
agreement is also very good, although our measurement is 
slightly larger than the prediction. For the Auckland clear filter. 



Fig. 5. First three basis functions of the PCA decomposition, 
computed for a set of atmosphere models that match the stellar 
parameters of OGLE 2004-BLG-482's source star 



parameter that defines the shape of the star's brightness profile 
is K = a2/0Ci. With our choice of PCA basis, all possible pro- 
files are obtained by varying K in the range [—0.162, 0.090], 
from the most peaked to the flattest limb-darkening profiles. 

We performed the OGLE 2004-BLG-482 analysis using the 
Heyrovsky ( 2003| l formalism for different combinations of data 
sets in a similar way as in Sect. |5] The results are presented in 
Fig. |6]for the combinations of data sets that were selected in 
the previous section (Fig. |4] panels with a black square in the 
upper right). 

As for the classical LLD law discussed in detail in the pre- 
vious section, we find a very good agreement between model 
predictions and our measurements. This shows that PCA LD 
provides an interesting alternative to analytical models of stel- 
lar brightness profiles. For applications where LD is not fitted 
(e.g.,I Kubas et al.|2008| l, it can be interesting to use PCA rather 
than LLD laws. On the other hand, the PCA LD law always de- 
pends on the set of selected model atmospheres. This could lead 
to discrepancies for instance if the parameter grid is too nar- 
row. In addition, any PCA LD law reflects the physics included 
in the construction of the particular atmosphere model (e.g., 
variants of ATLAS, MARCS, or PHOENIX models), which 
may not be ideally suited for the studied star In either of these 
cases, however, if the observational data are good enough, one 
may use the situation to one's own benefit. By trying different 
PCA LD laws and checking the quality of the fits and patterns 
in the residuals one can discriminate between different "can- 
didate" model atmospheres. To summarize, in limb-darkening 
modelling LLD has the advantage of simplicity and analyticity, 
while PCA LD has the advantage of providing better accuracy 
and flexibility. 



7. Conclusion 

We have analysed OGLE 2004-BLG-482, a relatively high- 
magnification single-lens microlensing event with notable 
extended-source effects, which was densely covered by our 
telescope networks. We derived precise Umb-darkening mea- 
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Fig. 6. PCA limb-darkening (PCA LD) coefficients K measured (crosses) and predicted (diamonds, Heyrovsky 20081 using the 
2-term PCA LD as explained in the text. Letters and colours have the same meaning as in Fig.|4] 



surements of the source star, a cool M giant, in particular in 
the / and R bands by combining the UTas and Danish data sets. 
No definitive conclusion could be made for Auckland data, af- 
fected by unknown systematics that prevented the data to be 
well-fitted along with other data; however, when the Auckland 
data are fitted alone, the derived limb-darkening agrees to 
model predictions, but with a large uncertainty. 

It provided us with the rather rare opportunity to di- 
rectly test model-atmosphere limb-darkening predictions for 
the source star This comparison was made possible because we 
could obtain high-resolution UVES spectra at VLT at a critical 
time thanks to the short activation of a ToO programme at VLT, 
from which we could precisely estimate the star's fundamental 
parameters. The source typing has been confirmed with good 
precision by our photometric diagnostic based on a calibrated 
colour-magnitude diagram of the field. We have performed a 
very detailed modelling to evaluate the impact of including data 
sets in the modelling process, and provide new diagnostics for 
future work. 

Very interestingly, the measured limb darkening agrees 
very well with model-atmosphere predictions, both when con- 
sidering linear limb-darkening profiles, or an alternative limb- 
darkening description based on a principal component analysis 
of ATLAS stellar atmosphere models. From this study, where 
the precision has been pushed to a high level, we conclude that 
this late M giant does not suffer from any clear discrepancy 
between limb-darkening model predictions and measurements, 
which has been pointed out for earlier K giants. Although it is 
based on the observation of a single event, it is very likely that 
the conclusion can be extended to similar late M giants. 
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